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Recent studies have been undertaken to investigate and develop aerodynamic models that 
predict aircraft response in nonlinear unsteady flight regimes for transport configurations. 
The models retain conventional static and rotary dynamic terms but replace conventional 
acceleration terms with more general indicial functions. In the Integrated Resilient Aircraft 
Controls project of the NASA Aviation Safety Program one aspect of the research is to apply 
these current developments to transport configurations to facilitate development of 
advanced controls technology. This paper describes initial application of a more general 
modeling methodology to the NASA Langley Generic Transport Model, a sub-scale flight 
test vehicle. 
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Nomenclature 

transfer function coefficients 

parameters of deficiency function 

normal-force slope of the body-wing and tail 

lift, normal force, and pitching-moment coefficients 

mean aerodynamic chord, ft 

frequency, Hz 

deficiency functions 

reduced frequency, n cf / V 

distance between aircraft c.g. and a.c. of the tail, ft 
pitch rate, rad/sec 

multiple correlation coefficient or coefficient of determination 

estimated standard error 

reference area, ft 1 2 

dimensional time constant, sec 

time, sec 

velocity, fps 

tail volume 

body-wing center-of-pressure position, fraction of c 

angle of attack, rad or deg 

mean (nominal) angle of attack, rad or deg 

amplitude of angle of attack, rad or deg 

downwash angle, rad or deg 

tail setting angle, rad or deg 

dummy integration variable 

non-dimensional time constant 

angular frequency, rad/sec 
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subscripts 

A = amplitude 

bw = body-wing 

t = tail 


superscripts 

B,T,W = body, tail, wing components 
BT, BW = body-tail, body-wing components 
BWT = body-wing-tail 



abbreviations 

c.g. = center of gravity 

a.c. = aerodynamic center 


I. Introduction 

F or conditions that include nonlinear or unsteady flight regimes, modeling and identification of transport aircraft 
in a classical wing-tail configuration brings new challenges to mathematical model postulation, wind tunnel and 
flight test experiment design, and aerodynamic parameter estimation. 

The problem of model formulation for a wing-tail combination 
in planar motion was addressed in Refs. 1-3. In all these cases only 
linear unsteady aerodynamics were considered and no application of 
the theory was made to experimental data. An extension to a 
nonlinear unsteady aerodynamic formulation and parameter 
estimation from wind tunnel data appeared in Ref. 4. In this analysis 
both steady and oscillatory data were used. The identified model 
was presented in several graphs comparing measured and estimated 
data. Unfortunately, no numerical values of estimated parameters 
and their accuracies were given. For that reason the same data were 
re-examined in Refs. 5 and 6, addressing mainly measured data 
inconsistency, low accuracy, non-physical values of parameter 
estimates, and model adequacy. 

In this paper the problem of system identification is considered 
again. This time, however, new sets of steady and oscillatory data 
for an airliner are used. These data were selected from extensive 
data sets produced at NASA Langley during the years from 2000 to 
2007. The main distinction between previous and recent sets of data 
is in the small amplitude oscillatory data. In recent tests, the number S = 5.90 ft 2 * , c = 0.915 ft, b = 8.85 ft 
and range of test frequencies were increased. The new data allow 

identification of a model with linear unsteady aerodynamics. In Figure 1. Model Geometry foi NASA 

future this model should form a basis for extended models including experimental sub-scale aircraft, 

nonlinear unsteady aerodynamics. 

II. Measured Data 

In order to obtain an extensive and detailed aerodynamic data base, three steady and dynamic low-speed wind 
tunnel tests were conducted during 2000-2001, 2003, and 2007, in the NASA Langley 14x22 Wind Tunnel, using a 
5.5% scale model representing a conventional twin-engine commercial airliner. This model has the same mold line 
specifications as the NASA Langley sub-scale flight test vehicle, referred to as the Generic Transport Model (GTM). 
The sketch and basic dimensions of GTM are given in Fig. 1. Tests were conducted mostly at a speed of 86 ft/sec, 
corresponding to a Reynolds Number of 0.5xl0 6 based on mean aerodynamic chord. The steady data were obtained 
for different configurations and an extensive range of angles of attack, sideslip angles, and control surface 
deflections. The dynamic tests included one degree-of-freedom forced oscillation around the pitch, roll, and yaw 
axes at different nominal values of the angle of attack, for different frequencies and amplitudes. Commanded 
oscillations were in the form of single sinusoids. An overview of the first two tests is presented in Refs. 7 and 8. 

For the following data analysis, only the normal force and pitching moment coefficient were considered. The 
steady data were taken from the 2003 experiment. They include the aerodynamic coefficients for body alone (B), 
body-wing (BW), body-tail (BT), and basic configuration (BWT) at angles of attack from -10° to 60° and for some 
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runs to 90°. Data for four configurations are plotted in Figs. 2 and 3. For the BWT case, measurements were made in 
both directions as the test rig moved up and down the angle-of-attack range. Consequently, at some angles of attack 
the plot will show two “+” symbols closely spaced. 

The oscillatory data were taken from the 2007 experiment for two configurations, BW and BWT, at nominal 
values of angle of attack in the range -10° to 50°, three amplitudes of 5°, 10°, and 15°, and typically six frequencies 
ranging from 0.01 to 0.86 Hz. During the test runs, data were obtained from 7 oscillation cycles for low frequency 
data, and 14 to 44 cycles for the remaining data. Oscillatory data were sampled at 250 Hz with 100 Hz analog filters. 
The resulting data were further filtered with a 4 Hz low pass digital filter to remove unwanted frequency content. 
The filter was run in both directions to ensure no phase error was added to the data. Initial cycles were discarded to 
avoid any startup transients in the data and five repeated cycles were used for data analysis to ameliorate the effects 
of measurement noise. 
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Figure 2. Normal-force coefficient for four 
model configurations. 



a, deg 

Figure 3. Pitching-moment coefficient for four 
model configurations. 


III. Analysis of Steady Data 

Selection of mathematical model structure for the analysis of steady data was influenced by the form of 
measured data (see Figs. 2 and 3) and by experience in modeling, as discussed in Refs. 4 to 6. In this study, model 
formulation assumes that the aerodynamics of the complete airplane can be formulated as a sum of two 
contributions, namely, the body-wing and the tail. This assumption is validated later in this section by showing that 

measured C BWT is equal to computed C BWT =C Bn + C J , where subscript a is A or m. The aerodynamic 
coefficients C N and C,„ can be expressed as 



C N {a) = C B N w {a) + S t C T N {a t ) 

(1) 


CJa) = C™' +C BW {a)x BW {a)-l,S t C T N {a r ) 

(2) 

where 

r T r BT r B 
L N ~ 

(3) 


a t = a + </> -£■(«) 

(4) 


Furthermore, S t = S t / S , £ t = £ 1 1 c , and x Bn is the distance between the body-wing aerodynamic center and the 
aircraft moment reference center as a fraction of mean aerodynamic chord. This distance can be computed from 
measured BW steady data and using Eq. (2) with the tail terms removed. C N or C m without superscript refers to 

the complete BWT configuration, C BWT or C BWT , respectively. The downwash angle at the tail, s(a ) , follows 
from solving the following two equations for a t at selected a. 
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( 5 ) 


Cl(a) = C» T {a)-C»{a) 

Cl(a) = -l t S t C T N (a t ) (6) 


The variables x m ( a ) and e(a) are plotted in Figs. 4 and 5. 
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Figure 4. Distance from moment reference 
center to aerodynamic center. 



Figure 5. Downwash angle at the tail. 


In order to allow steady data computations at any angle of attack, the measured data were fitted by high-order 
polynomials in a. Data compatibility was checked by comparing measured and computed C N and C„, using Eqs. (1) 
to (4). The results are presented in Fig. 6. The measured and computed coefficients are in very good agreement thus 
substantiating the use of superposition in forming model equations. The additional check confirmed the equality in 
Eq. (6) and the validity of £ t = 4.02. This value is obtained from the model geometry. 

Force contributions from the isolated wing and isolated tail can be computed from steady data using Eqs. (3-4) 
and from the expression 


r W _ /-tBWT r BT 
L N ~ c N L N 


( 7 ) 


and similarly for moment contributions, and . Steady data analysis was completed by considering C L (a) , 


plotted in Fig. 7, since this coefficient shows pre- and post-stall regions more clearly than C N (a) . 





a, deg 

Figure 6. Measured and computed coefficients 
for complete aircraft. 


a, deg 

Figure 7. Computed lift coefficient for complete 
aircraft from normal force and axial force. 


IV. Model Formulation with Linear Unsteady Aerodynamics 

Development of linear model equations for an aircraft in one planar motion follows from the results in Refs. 3, 9, 
and 10. For this study the normal force and pitching moment can be considered as functions of angle of attack and 
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pitch rate, i.e., C a (t) = C a (a(t),q(t )) , where a = N or m. Then the aerodynamic model equations can be formulated 
as 


c a (0 = c a (0) + } C aa (t - z)d{z)dz + 1 C aq (t - z)q(z)dz (8) 

o 2K 0 

where C a (0) is the value of the coefficient at initial steady-state conditions, and C aa (t ) and C a ^ ( t) are the indicial 

functions defining the responses in C a to unit step in a and q, respectively. Replacing the second indicial function 
in Eq. (8) by its steady value C a (oo) , Eq. (8) is simplified as 


C a (t) = C a ( 0 ) + — Cr (°°)<?(0 + 1 C aa ( t - z)d(z)dz 
IV H o 


(9) 


When the deficiency function 


F aa {t) = C aa ^)-C aa {t) 


(10) 


is substituted into Eq. (9), the equations for the aerodynamic coefficients take the form 

C N (t) = C, v (0) + C N (oo )a(t) + -^—C N (oo)q(t) - 1 F n (t - z)d(z)dz (11) 

Z.V * Q 

Cm (0 = C m (0) + C,„ a (oo )a(t) + ^C mq (°°)?(0 -\ F m a it - r)d(z)dz (12) 

where C Na (oo) , C ma (oo) , C N (oo) , and C m (oo) are the partial derivatives with respect to a and q, respectively, 

of the coefficients C N and C m in steady flow conditions. 

The indicial functions in Eqs. (11) and (12) include combined responses of the wing, tail, and interference effects 
between the lifting surfaces. With the assumption of linear aerodynamics the resulting indicial functions are given as 
a sum of four components. For the indicial function, C Na (/ ) , these components are as follows: 

1 . the response in the normal force of the wing to a unit step in angle of attack of the wing, while the angle 
of attack of the tail remains at zero; 

2. the response in the normal force of the tail to a unit step in angle of attack of the tail, while the angle of 
attack of the wing remains at zero; 

3. the response in the normal force of the tail to a unit step in angle of attack of the wing, while the angle 
of attack of the tail remains at zero; 

4. response in the normal force of the wing to a unit step in angle of attack of the tail, while the angle of 
attack of the wing remains at zero. 

The first two components represent the response of an isolated wing and isolated tail, respectively. The third 
component expresses the lift on the tail due to a change in the downwash induced by lift of the wing. It is usually a 
significant contribution to the resulting pitching moment of an aircraft. For the following data analysis the second 
and fourth components of the indicial function were assumed to be small relative to total forces or moments and 
therefore were neglected. 

Contributions of the body-wing and tail to coefficients C N and C m are described in the relations 

C N =C B N W +S t C T N (13) 
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( 14 ) 


r - r , BW n BW j 7r T 

'-'m ^ m a.c. X '-'iV 


where C m ^ < is the pitching moment about the aerodynamic center of the body-wing combination, x m is the 

distance between the BW aerodynamic center and the aircraft center of gravity or moment reference center, and tail 
volume, V t = £ t S , . Combining Eqs. (11-14) yields 


C N U) = C N ( 0) + C Na (oo)a(O + y - C Nq (co) q (t) - j [f™ (t- r) + S,F^ (t - r)]a(r)rfr (15) 

AV 0 

Cm (0 = C ma c + C ma (00 )a{t) + A-C mq ( 00 MO - \[x? w F™ (t- r) + ^ (# - r)]«(r )rfr (16) 

where the deficiency functions are defined as 

T?BW u\ /^BW / \ /^BW />\ 

F ,v 0 (0 = Qv a (<*>)- C Na (t) (17) 

F N a (!) = C N a <*>)- C N a (t) ( 18 ) 


The steady terms in Eq. (15-16) include the aircraft stability derivatives, , C„ a , C N , and C m . Indicial 
flinctions are reduced to these terms under steady flow conditions. The first two derivatives are usually formed as 

C N a =a bw + a t S t (l-ds/da) (19) 

C ma = x BW a bw -a t V t (l-d £ /da) (20) 

In these equations a bw = 8C BW / da and a, = dC ^ / da, are the normal force slopes of the body-wing and tail, 
respectively. The derivatives C Nq and C mq represent aerodynamic effects that accompany rotation of the aircraft 

about its lateral axis while a remains constant. Both body-wing and tail are affected by the rotation but the body- 
wing contribution is usually negligible. Thus considering a contribution of the tail only, the damping parameters are 
approximated as 

C Nq =2a t V t (21) 

C mq =-2a t l t V t (22) 

For parameter estimation from measured data the form of the deficiency functions in Eqs. (15-16) must be 
specified. Experience with identification of tailless aircraft indicates that a simple exponential function 

Fjf a (f) = ae~ br (23) 

T . 

can be a good approximation. In Ref. 3 the model for indicial function C N ( t ) is based on linear aerodynamics and 
a combination of three effects: (1) change in C^ n due to change in a, (2) change in fat the tail due to a change in 
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Ctf W and (3) by a change in C] due to change in s. This model, however, is not suitable for parameter estimation 
due to its complexity. 

T 

A simplified form of C N was proposed in Ref. 4 assuming that the propagation of the downwash is governed 
by a first order differential equation 


Te = -e + s( oo ) (24) 

leading to an exponential form of the deficiency function. In this paper a simple form was used for the deficiency 
function 


Fl{t) = 


[constant for 0 < t < At ) 


lo 


for t elsewhere] 


(25) 


This formulation is based on the Cowley-Glauert explanation of the downwash delay . They assumed that a sudden 
change in the angle of attack of the wing results in a sudden change in the downwash. Downwash from the wing is 
delayed from reaching the tail by the time A t = l t IV . Then the contribution of the tail to the normal force 1 1 is 

a, Qs 

A C N — -a,S t As = a t S t — — a 
V da 

or further simplifying as 


A C J N =2 a t £ t S t 


? ds dc t dc 


da 2V 


= C 


Na 


2V 


(26) 


and to the pitching moment 


dC T =2at SV— ^=C T — 
m 1 ' ‘ 1 da 2V 2V 


(27) 


These formulations replace the second integral in Eqs. (15-16) by a constant. From Eqs. (15-16) and Eqs. (26-27) 
the model equations for the increments in the normal force and the pitching moment take the form 


C N (t) = C N (0) + C N a(t) + ^-C' N q(t)-a\e h(t T) d(r)dr 

IV y Q 


(28) 


CJt) = C_ +C m a(t) + ^C'q(t)-a\x BW e- b ^d(T)dT 


2V 1 


(29) 


where 


C Nq =C Nq+ C' N . 


(30) 


c =c 

tiiq n 


+ct 


(31) 


In model Eqs. (28-29) there are six parameters. Four parameters are the aerodynamic derivatives, C Na , C,„ a , C' N , 
and C' , and two parameters, a and b, characterize the deficiency function. 
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By introducing 


i) = je h{t r ^d(r)dr 
o 

and omitting the bias term, the state-space format of Eq. (28) is 

f](t) = -bij(t) + d{t) 


( 32 ) 


(33) 


C N (t) = C Na cc(t) + — C' Nq q(t) ~ arj(t) 


(34) 


Applying the Laplace transform to Eq. (33-34), the transfer function for the normal-force coefficient is obtained as 


Cjy (s) As + Bs + C 


a(s) 


■ +b 


(35) 


where s is the Laplace transform variable and 


A = — Ca 


2V 


N„ 


B = C N -a + b — C' N 
"a 2V q 

C = bC Nn 


(36) 


are the transfer function coefficients. 


V. Model Identification 

As mentioned previously, measured oscillatory data were obtained at different angles of attack, amplitudes, and 
frequencies for two configurations: the basic model (BWT configuration) and the model without tail (BW 
configuration). In this study data analysis was primarily focused on data with a A =10° because these data formed a 
set with twenty nominal values of a, from -10° to 50°, and seven frequencies, from 0.05 to 0.86 Hz. For the set of 
data with a A = 5°, the number of frequencies was reduced to six and the number of a 0 to twelve. Data with a , = 15° 
were not analyzed for this paper but will be used for future research in nonlinear modeling. Model identification for 
this study includes an assessment of model linearity and its dependence on angle of attack, amplitude, and 
frequency. Model identification is completed by parameter estimation and model validation. 

A. Harmonic Analysis 

A method of harmonic analysis 12 was applied to measured aerodynamic force and moment coefficients to allow 
estimation of in-phase and out-of-phase coefficients. Five cycles of forced-oscillation data were analyzed for each 
frequency. A mathematical model was postulated as 


m m 

C a (t) = Aq + X A - cost jcat)+ X B: sin(y'ot) a = N ov m (37) 

j = 1 7=1 

where co = 2nf and A 0 , A ,, B lt A„„ B m are the Fourier coefficients. The analysis provided estimates of these 
coefficients, their standard errors, and the coefficient of determination, R . This coefficient indicates how much 
information in the data is explained by the model 1 '. For a model with linear aerodynamics, the in-phase and out-of- 
phase components of C a (t ) , i.e., C a ^ and C aq , can be expressed in terms of coefficients A , and B I as 
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a = N or m 


(38) 


C a = — =C a -k 2 C a . 

“a aA a a “q 


C a = — l -=C a +c a . 
‘I ka A <1 a 


a = N or in 


(39) 


where k = ncf IV is the reduced frequency. Also shown in Eqs. (38-39), are the commonly expressed relationships 
of these coefficients to the steady flow damping terms and acceleration terms. 




a, deg 

Figure 8. Normal-force coefficient from static 
and oscillatory data,/ = 0.08 Hz, a 4 = 10°. 


a, deg 

Figure 9. Normal-force coefficient from static 
and oscillatory data,/ = 0.86 Hz, a A = 10°. 


Examples of the measured normal-force coefficient from steady and oscillatory tests are shown in Fig. 8 and 9. 
Five cycles of oscillatory data are presented for a 0 at -5°, 14°, 30°, and 45°, and for /= 0.08 and 0.86 Hz. Fig. 9 
highlights the effect of increased steady-flow damping at higher / and do > 10°. The departure of the oscillatory data 
from a regular ellipse indicates a departure from linear aerodynamics. 

The estimated parameters C aa and C a ^ for a = N or m, are presented in Fig. 10 and 11 for BWT configurations 

and a 4 = 10°. The estimates of C Na and C„, a change very little with frequency and their values agree closely with 

derivatives C Na and C m ^ estimated from steady measurements. On the other hand, parameters C N ^ and C m ^ 

exhibited strong frequency dependence for / < 0.43 Hz and 10° < ag< 50°. Both figures also include estimates of R 2 . 
These values indicate that a linear aerodynamic model may be adequate for -10° < «» < 30°. Figure 12 includes only 
the out-of-phase components of C N and C„, for BW configuration and a , = 10°. Figs. 10-12, show the contribution 
of the tail to aircraft damping. 

In order to explain the variation of the in-phase and out-of-phase components with frequency these components 
were expressed in terms of deficiency functions 14 ’ 15 , as in Eq. (23). In Refs. 14-15, the model for these components 
has the form 


r 2 k 2 


C “«- C “«- a X + z 2 k 2 


(40) 


c„ =C„ -a- 


1 + T 2 k 2 


( 41 ) 
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where subscript a = N or m. The second terms on the right side of Eqs. (40-41) can be considered as an unsteady 
counterpart of k 2 C a and C a& . For the model Eqs. (40-41), there are four unknown parameters ( C a ^ , C a ^ , a , T\ ) 


which can be estimated from C n and C„ . 

c, a “ q 



Figure 10. Harmonic analysis for normal-force Figure 11. Harmonic analysis for pitch-moment 

coefficient, BWT configuration, ^=10°. coefficient, BWT configuration, a^=10°. 


To demonstrate the effect of amplitude on the normal-force coefficient, a plot of five cycles of C N (t)v s a is 
provided in Fig. 13 for a A = 5°,/= 0.86 Hz, and a n at 0°, 10°, 20°, and 30°. As in the a A = 10° case, the steady 
values of C N were included for comparison. The small amplitude oscillatory data forms shapes much closer to a 
regular ellipse, except near the stall break. This supports an assumption of linear aerodynamic responses for these 
conditions. 




Figure 12. Out-of-phase component of normal 
force and pitching moment, BW, a A = 10°. 


a, deg 

Figure 13. Normal-force coefficient from static 
and oscillatory data,/= 0.86 Hz, BW, a A = 5°. 


Figs. 14-15 contain out-of-phase components, C N ^ and C m ^ , 


for both BWT and BW configurations with a A = 
5°. Some differences are visible when comparing these results with the estimates at a A = 10°. The main difference is 
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the effect of frequency on C m for 10° < a 0 < 30°. Explanation of this effect, amplitude effects, and other 
differences will be covered in future research. 



Figure 14. Out-of-phase components for normal- Figure 15. Out-of-phase components for normal- 
force and pitching moment, BWT, a A =5°. force and Ditching moement. BW. a A = 5°. 


If the principle of superposition applies to forced-oscillation data, then the damping parameter of the tail alone 
can be computed as 


t; t _ t^bwt t;bw 

L N q ~^N q <~N q 


(42) 


Parameter C)y is plotted in Fig. 16 against frequency 

and four values of cto- This plot indicates little 
frequency effect for / > 0.2 Hz, however, for lower 
frequencies the value is rapidly increased. Similar 

behavior was observed for parameter C m ^ . These 
changes and the effect of introducing deficiency 

T 

function F N are not explicitly included in the 

current model formulation, see Eqs. (28-29), and have 
not been studied. 


B. Parameter Estimation 

For estimation of aerodynamic derivatives and 
indicial (or deficiency) function parameters, two 
approaches are used: 1) a two step linear regression of 
Ref. 12 or 2) a variation of an output error method 
explained in Ref. 13 and implemented in Ref. 16. The 
two-step linear regression is based on model Eqs. (40) and (41). This method estimates r> in the first step as the only 
unknown parameter. Then, in the second step, zq is treated as a known parameter and the remaining three 



0.4 
f, Hz 

Figure 16. Out-of-phase components for tail 
alone, a A = 10°. 


parameters, C Na , C N and a , are estimated. This method works very well provided that a distinct linear 
relationship exists between C Na components and C N for different frequencies at a given a n . This is not the case for 
the present experiment as follows from the lack of frequency dependence in C a ^ , shown Figs. 10 and 11. In this 
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study an output error based method 9 was used for estimation of parameters in state-space model Eqs. (33) and (34) 


to utilize harmonic data. In this case, inputs, a(t ) , 
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Figure 17. Variation of unsteady term parameters, t and a, with frequency. 


that allowed use of the time-domain output error algorithm 
d(l) and output 
measurements C a (i) 
were stacked for six 
frequencies to capture 
as much information 
and frequency content 
from the data as 
possible for each 
selected cto- The a(t) 
input time history was 
obtained by 

differentiating and 
smoothing five cycles 
of measured a(i) . 

Parameter estimation 
results were obtained 

at four values of a (h two configurations, BWT and BW, and two amplitudes, a A = 5° and 10°. The results show that 
estimates of Cy^ were very close to C N from harmonic analysis and C N from steady measurements. The average 

standard error of estimated parameters C Na was 0.0025. Estimates of C v ^ with an average standard error 0.68 came 

close to C N for the highest frequency covered by the 

experiment. Estimates of parameter b were converted 
into estimates of the time constant, T = l/b . This 
parameter is plotted in Fig. 17 against a for two 
configurations and two amplitudes. The results indicate a 
gradual increase of time constant with a. For better 
clarity, parameter error bars were not included in Fig. 

17. Variation of estimated parameter a for different 
configurations and angles of attack is also shown in Fig. 

17. Trend in the data follows an expected variation from 
zero (in low a regions without unsteady effect) and back 
toward zero for higher a where the steady model is 
adequate. 

Representative contributions from each term in 
the C N (t) output equation are shown in Fig. 18 for a 
selected frequency. The main contributors are the static 
and steady flow damping term. The unsteady term 
provides a relatively small component to the total 
response for each frequency and each of the four a 0 cases considered. 



-10 


0 


10 


30 


40 


50 


20 
a, deg 

Figure 18. Contribution of three terms to 
normal-force output equation, BWT, cca=10°. 


C. Model Validation 

For validation of the model as a predictor, a comparison of measured and predicted data was used. Fig. 19 show 
the variation of Cw with a, at /= 0.86 Hz, two nominal values of cto (14°, 26°) and a A = 10°. In both cases the shape 
of the predicted data forms a regular ellipse reflecting the underlying linear aerodynamic model structure. The 
measured C v shows some deviation from a regular ellipse, or linear behavior, in the pre-stall and stall region. At a ( , 
= 26°, the behavior of the data supports using a linear model. Considering the data for a A = 5°, the predicted data are 
very close to the measured data for both a,, conditions. However, the shape of measured data for a n = 14 deg is still 
somewhat affected by the aerodynamics of the stall region. 
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Figure 19. Comparison of measured and predicted Figure 20. Comparison of measured and predicted 
normal-force coefficient, BWT, f = 0.86 Fiz, au=10°. normal-force coefficient, BWT, f = 0.86 Hz, a A = 5°. 

VI. Concluding Remarks 

This paper presents aircraft model identification from wind tunnel data of a transport aircraft and its components: 
body-wing-tail, body-wing, and tail. The measured data include results from steady and oscillatory tests over an 
extended range of angle of attack, frequency, and two amplitudes. 

Using the steady data it was shown that the measured normal force and pitching moment for a complete model 
are identical to those obtained as a sum of body-wing and tail contributions. Accepting a superposition principle, it 
was possible to estimate the distance from a moment reference center and the body-wing aerodynamic center, as 
well as the downwash angle at the tail. It was also possible to prove a consistency between normal-force and 
pitching-moment data, and to compute the normal force of isolated wing and isolated tail. 

Static stability derivatives, C aa , estimated from oscillatory data were consistent with calculated values using 

high-order polynomials in a fit to static measurements. 

The analysis of oscillatory data was preceded by formulation of a mathematical model for aircraft in planar 
motion with one degree of freedom and with linear unsteady aerodynamics. This model, in the state space form, was 
later used for parameter estimation. The data for this approach were obtained from measurements with two 
configurations (complete model and model without tail), seven frequencies, and two amplitudes of 5° and 10°. 

For harmonic analysis of the data, a mathematical model was defined by equations of motion representing a 
steady harmonic solution. Model parameters included the in-phase and out-of phase components of normal force and 
pitching moment and two parameters of an indicial function: time constant, T, and a gain, a, characterizing the level 
of contribution from the unsteady aerodynamic terms. Estimates were completed by including parameter standard 
errors and the coefficient of detennination for normal force and pitching moment models. 

Harmonic analysis indicated that a linear unsteady model is adequate for angles of attack from -10° to 30° for 
both amplitudes, although some limited nonlinearities present in the pre- and post-stall regions slightly reduced the 
effectiveness of the linear model. The in-phase components are very close to results from steady data and the out-of- 
phase components at high frequencies are close to steady damping term. The indicial-term time constant is 
increasing with an increase of angle of attack. The second parameter of the indicial function correctly indicated the 
region where the unsteady effect is close to zero. 

The resulting model was validated by its prediction capabilities. Measured and predicted normal-force 
coefficients at different nominal values of the angle of attack and amplitudes were compared. As expected the linear 
unsteady model with parameters estimated from low amplitude data was a better predictor than the model with 
parameters from experiment with higher amplitude. 

Future work will expand model identification of the NASA Generic Transport Model for lateral-directional axes 
and investigate higher amplitude and higher a regions where the nonlinear unsteady model is required. 
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